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ABSTRACT 

In this paper, we report the detections of stellar variabilities from the first 2-year observations of sky 
area of about 1300 square degrees from the Tsinghua University-NAOC Transient Survey (TNTS). A 
total of 1237 variable stars (including 299 new ones) were detected with brightness < 18.0 mag and 
magnitude variation > 0.1 mag on a timescale from a few hours to few hundred days. Among such 
detections, we tentatively identified 661 RR Lyrae stars, 431 binaries, 72 Semiregular pulsators, 29 
Mira stars, 11 slow irregular variables, 11 RS Canurn Venaticorum stars, 7 Gamma Doradus stars, 5 
long period variables, 3 W Virginis stars, 3 Delta Scuti stars, 2 Anomalous Cepheids, 1 Cepheid, and 
1 nove-like star based on their time-series variability index J s and their phased diagrams. Moreover, 
we found that 14 RR Lyrae stars show the Blazhko effect and 67 contact eclipsing binaries exhibit the 
O’Connell effect. Since the period and amplitude of light variations of RR Lyrae variables depend on 
their chemical compositions, their photometric observations can be used to investigate distribution of 
metallicity along the direction perpendicular to the Galactic disk. We find that the metallicity of RR 
Lyrae stars shows large scatter at regions closer to the Galactic plane (e.g., —3.0 < [Fe/H] < 0) but 
tends to converge at [Fe/H]~ —1.7 at larger Galactic latitudes. This variation may be related to that 
the RRAB Lyrae stars in the Galactic halo come from globular clusters with different metallicity and 
vertical distances, i.e. Ool and OoII populations, favoring for the dual-halo model. 

Subject headings: stars: variables: general; stars: variables: RR Lyrae; binaries: eclipsing; 
stars:evolution 


1. INTRODUCTION 

Time-Domain Astronomy (TDA) has emerged as a key 
field of current astronomy and astrophysics (Grindlay et 
al. 2012). This offers new routes to astrophysical under¬ 
standing of stellar and galactic activities such as variable 
stars, novae, supernovae, gamma-ray bursts, AGN and 
quasars etc.. The Tsinghua University-NAOC Transient 
Survey (TNTS) is a four-year project designed to search 
such transients (Zhang et al. 2015). The observations 
were carried out in a clear filter with a 0.6-rn Schmidt 
telescope located at Xinglong station of NAOC. Under a 
typical weather condition at Xinglong station (e.g., with 
a seeing of about 2 arcsec), this telescope and the CCD 
system can reach a detection limit of about 19.5 mag (3cr) 
in the clear filter for an exposure of 60s. The relatively 
short cadence (e.g., 3-4 days) and broad synoptic cover¬ 
age of the sky areas provide an excellent chance to search 
for and study variable stars of different time scales. In 
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this paper, we present the result of searching for variable 
stars from these survey data. 

This work presents the result of searching for variable 
stars from the TNTS data accumulated during the period 
from Oct. 2012 to Aug. 2014. The observations and data 
reductions are described in Section 2. The time-series 
photometry from TNTS images is presented in Section 
3. The classifications of variable stars and analysis are 
given in Section 4. We summarize in Section 5. 

2. OBSERVATIONS AND DATA REDUCTIONS 

The TNTS is conducted with a 0.6-m Schmidt tele¬ 
scope equipped with a 4kx4k CCD. The CCD has a pixel 
scale of 1.3 arcsec pixel” 1 and a field of view (FoV) of 1.5 
degrees x 1.5 degrees. It’s designed to operate for four 
years starting from October 2012. This survey covers a 
sky area of ~1800 square degrees with Galactic latitude 
|6| > 10° and 0° < S < 60°. It usually takes about 2 
minutes to get an image for a specific sky zone, includ¬ 
ing 60-s exposure for the selected target, 22-s readout 
time, and 30-s movement and stabilization for the tele¬ 
scope. In order to efficiently rule out the cosmic rays and 
moving objects, we took two exposures for each sky zone 
with a temporal interval of 1-1.5 hours and the whole 
survey area can be visited every 3 to 4 days (see details 
in Zhang et al. 2015). 

About 85,000 images were accumulated from our sur¬ 
vey performed during the period from Oct. 2012 to Aug. 
2014, which cover about 1800 square degrees of the north¬ 
ern sky. In this study, however, we only chose the fields 
with repeated observations > 40 times to detect variable 
stars. This selection criterion results in that the observa¬ 
tions of a sky area of ~ 1300 square degrees can be used 


2 


Yao et al. 



Fig. 1.— The photometric error versus magnitude in a clear band 
for a typical image of the field SE65 from the TNTS. 


for such an analysis. The TNTS images were reduced us¬ 
ing the IRAF 7 standard procedure, including the correc¬ 
tions for bias, flat field, and removal of cosmic rays. The 
positions of those images were then registered to each 
other using the Guide Star Catalogue (GSC) (Jenkner 
et al. 1990), with an accuracy of 0.5 arcsec (Zhou et al. 
2003). The SExtractor (Bertin & Arnouts 1996) was used 
to perform aperture photometry on the chosen images. 
The radius of aperture photometry was set as 5 pixels 
(6.8 arcsec), with a sky annulus size of 24 pixels (31.2 
arcsec). 

Since the stars on the TNTS images are not very 
crowded, we selected the parameter ’’BEST OF MAG” 
of SExtractor as output results. The detection thresh¬ 
old was set as 2<r above the sky background. To alle¬ 
viate the effect of nonuniform background on the pho¬ 
tometry, i.e., the moonlight and clouds etc., we set the 
BACKGROUND TYPE as LOCAL. A range of 5000- 
17000 stars can be detected on the TNTS images. In our 
study we take the typical sky field SE65 as an example 
for illustration of photometric error. For stars brighter 
than 17.0 mag, the photometric error is generally less 
than 0.02 mag according to the output result from the 
SExtractor (see Fig. 1). 

3. TIME-SERIES PHOTOMETRY 

To get the time-series photometry of each star, we feed 
the matching tool SRPMatch with photometry catalogs 
produced by SExtractor for two images of the sky field. 
SRPMatch is a command line tool embedded in the Swift 
Reduction Package (hereafter SRP 8 ), which can be used 
to match photometry catalogues. We set the match tol¬ 
erance as 1.5 arcsec (~ 1 pixel). For each field, one is the 
reference image with the maximum number of detected 
stars and the other is the random image from the same 
field. We calibrated our time-series photometry by the 
following two steps: magnitude calibration and rescaling 
of the photometric errors. 

' IRAF, the Image Reduction and Analysis Facility, is dis¬ 
tributed by the National Optical Astronomy Observatory, which is 
operated by the Association of Universities for Research in Astron- 
omy(AURA), Inc. under cooperative agreement with the National 
Science Foundation(NSF). http://iraf.noao.edu/ 

8 https://pypi.python.org/pypi/SRPAstro 


3.1. Magnitude Calibration 

To calibrate magnitudes of images taken at different 
epochs, it is necessary to calculate the magnitude offset 
of each image relative to the reference image. Fig. 2 
shows the magnitude offsets measured for the sky field 
SE65 (which is the same as that used in Fig. 1) and its 
reference frame. One can see that the photometric error 
progressively increases for stars fainter than 14.0 mag. 
On the other hand, the stars brighter than 12.5 mag tend 
to be at the edge of saturation on the TNTS images. We 
thus chose those stars with magnitudes between 12.5 mag 
and 14.0 mag to determine the magnitude offset for each 
image. To mitigate the effect of nonuniform background 
due to moonlight or obscuration by clouds, we divided 
every image into 16 small subsections and then calculated 
the offset of each subsection to perform the photometric 
calibrations. Outlier rejections of 3c were applied itera¬ 
tively to make sure that no variables were involved in the 
calculations. The median value of the magnitude offsets 
of 16 subsections is almost the same as that of the whole 
image, e.g., A M zp = 0.010 ±0.001 for the sky field SE65. 
The uncertainty in the magnitude offset is so small (rela¬ 
tive to the photometric errors from SExtractor) that we 
neglected it in the following calculations. 
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Fig. 2. — Magnitude offsets between the same image of the field 
SE65 used in Fig. 1 and the reference frame of the same field. 


Previous studies indicate that unfiltered magnitude is 
approximately equivalent to that in the R band of PPMX 
system (Position and Proper Motions extended, Roser et 
al. 2008), with an uncertainty of 0.2-0.3 mag, depending 
on the number of flux-calibrated stars available in the 
field (e.g., Li et al. 2003, 2011). For each field of the 
TNTS, we match the reference image with the PPMX 
catalogue to choose appropriate stars to calibrate the 
unfiltered magnitudes. Owing to the large FoV of the 
TNTS images, the number of stars that can be used for 
the flux calibration is relatively large. The median un¬ 
certainty in the photometric calibration for the 600 fields 
used in our study is ~0.010±0.005 mag, which is much 
smaller than that from Li et al. (2003). 

3.2. Rescaling the Measured Photometric Errors 

For some reason, the photometric errors given by SEx¬ 
tractor might be overestimated or underestimated. We 
would either miss real variables, or select bogus vari- 
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ables as real ones (Kaluzny et al. 1998). Assuming that 
the photometric errors roughly follow the Gaussian dis¬ 
tribution, the logarithm of the x 2 /Ndof of stars with 
constant luminosities should be close to zero (see the 
right panel of Fig. 6 from Kaluzny et al. 1998). The 
X 2 /N dof is defined as 


X 


2 


^ in — 


( 1 ) 


N DO f =N-M. (2) 

Where m is the weighted mean value (the weighting fac¬ 
tor is the inverse square of the error a m ) of the time- 
series magnitude m of a star, N is the number of mea¬ 
surements, and M is the number of parameters (here 
we let it be one). Take the field SE65 as an example. 
Fig. 3 shows the plot of the logarithm of the y 2 /Ndof 
against weighted mean magnitude for all of the stars. 
One can see that the brightest stars (~ 12.0 mag) have 
X 2 /N D of ~ 100, and their photometric errors are thus 
seriously underestimated. The possible reasons for this 
discrepancy includes underestimation of the flat fielding 
error and imperfect aperture photometry etc. To elim¬ 
inate this inconsistency, the magnitude error a is mul¬ 
tiplied by an appropriate scaling factor F, following the 
method proposed by Kaluzny et al. (1998). 
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Fig. 3.— The Log(x 2 /Ndof) distribution as a function of the 
weighted mean magnitude for stars in the field SE65 before rescal¬ 
ing their photometric errors. 


In detail, we first binned the data X 2 /Ndof at an in¬ 
terval of 0.1 magnitude. In the logarithm plot, we used 
a second-order polynomial to fit the data for stars 
with brightness < 13.5 mag. For stars with brightness 
between 13.5 mag and 18.0 mag, we used a linear function 
to fit. While those stars with weighted mean magnitudes 
fainter than 18.0 mag were discarded in the fit due to rel¬ 
atively larger photometric errors. We then applied the 
scaling factor F to rescale the photometric error of each 
star. The distribution of the rescaled photometric errors 
for stars in the field SE65 is shown in Fig. 4. 

4. VARIABLES FROM THE TNTS 
4.1. Searching for Variable Stars 
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Fig. 4.— The same diagram as the Fig. 3 but after rescaling the 
photometric errors. 


The variable candidates were initially selected using 
Welch-Stetson variability index J s (Stetson 1996). We 
used VARTOOLS (Hartman et al. 2008) to compute the 
J s index based on the rescaled errors of photometry. For 
each star one can calculate the variability index J as 


E n -> 

k =1 W k 


( 3 ) 


where sgn() is the sign function. We considered n pairs 
of observations, and each has a weight uy , 


/ if i(k) ± j(k), 

{ 6 i(k) - f> if *( fc ) = 3 


( 4 ) 


where Pk is the product of the normalized residuals of 
the two paired observations i and j , and 


c / n m — m 

6= \ -7- ( 5 ) 

V n - 1 a m 

represents the magnitude residual from the mean value 
scaled by the standard error (Stetson 1996). In detail, n 
is the total number of observations contributing to the 
mean. The final variability index is 

Js = ( 6 ) 

'Wmax 

In our calculations, we considered the two observations 
taken on the same night as a pair and give each of them 
a weight of 0.5. Otherwise a weight of 1.0 was given if 
there is only a single observation. 

As an example, we show in Fig. 5 the distribution of 
J s obtained for stars in the field SE65. As it is expected 
most stars have J s values close to zero. The mean J s 
value of 600 sky fields involved in this study is 0.07. For 
each field, we chose those stars with J s larger than the 
mean value by 3<r as candidates of variable stars. 

We searched for periodic variabilities among candi¬ 
dates using Lomb-Scargle (L-S) method (Lomb 1976; 
Scargle 1982; Press & Rybicki 1989; Press et al. 1992) and 
Phase Dispersion Minimization (PDM) method (Stelling- 
werf 1978). The L-S method applies the statistical prop¬ 
erties of least-square frequency analysis of unequally 
spaced data on a series of test periods. We set test pe¬ 
riods ranging from O.ld to 600d with a bin size of O.ld. 
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Fig. 5.— Distribution of the J a index calculated for stars in the 
sky field SE65. 

We determined the most likely period for each star in 
terms of the power value given by the L-S method. Then 
we performed PDM method to double check the periods 
detected from the L-S method since the PDM method is 
useful for data sets with gaps, non-sinusoidal variations, 
poor time coverage or other problems that would make 
Fourier technique unusable. The folded light curves were 
then established using the period from the L-S/PDM or 
twice of that for a binary if necessary. The period derived 
from the PDM method was only adopted if the phased 
light curve folded by the PDM has smaller scatter than 
that from the LS method. We examined the light curve 
and phase diagram for each candidate to pick out vari¬ 
ables and throw out those showing alias variations (e.g., 
one day, half day, one third day and so on) and having 
large dispersions. Finally we found 1237 variables with 
magnitude variations >0.1 mag and brightness < 18.0 
mag based on the L-S method. 

All of these candidates had been cross-checked with 
those listed in the VSX database (the International Vari¬ 
able Star Index, Watson et al. 2014), which is created by 
amateur astronomers of American Association of Vari¬ 
able Star Observers (AAVSO) 9 . We found that 938 out 
of these 1237 objects have been included in this database 
with 438 ones from the Catalina Sky Survey (CSS) sam¬ 
ple (Drake et al. 2014) and 106 ones from the Lincoln 
Near-Earth Asteroid Research (LINEAR) survey sample 
(Palaversa et al. 2013), and the rest 299 ones can be con¬ 
sidered as new variables. Fig. 6 shows the distribution 
of these variables, in which they are displayed accord¬ 
ing to the Galactic coordinates. The parameters of the 
first 20 ones of variable stars from variable star candi¬ 
dates table are listed in Table 1. Column 1 and 2 list 
the survey field and the ID number of the field; Column 
3 and 4 give the right ascension and declination of the 
relevant field; Column 5 gives the mean unfiltered magni¬ 
tude (calibrated with the R-band magnitudes of PPMX 
catalogue); Column 6 gives the line-of-sight extinction 
based on Schlafly & Finkbeiner(2011); Column 7 gives 
the amplitude of variation (magnitude from maximum 
to minimum); Column 8 gives the variability index; Col¬ 
umn 9 and 10 give the periods of the variables from our 
measurements and the corresponding method; Column 

9 AAVSO: http://www.aavso.org/vsx/ 


11 gives the period from the VSX; Column 12 gives the 
minimum light of the variables detected in the observa¬ 
tion season; Column 13 gives a tentative classification of 
the variable star whenever possible; Column 14 shows 
the additional information about the variable stars, in¬ 
cluding previous identifications from the VSX database 
(e.g., CSS or LINEAR). Note that, for those known vari¬ 
ables, the periods estimated from the TNTS are highly 
consistent with those given by the VSX database except 
about thirty ones (see Figure 7). The phased light curves 
of new periodic variable stars are shown in the appendix. 
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Fig. 6.— Distribution of the TNTS variable stars in the Galactic 
coordinates. 1 0 



Fig. 7. — The measured periods from the TNTS data versus those 
from the AAVSO VSX for the variables in common. 


4.2. Classification of Variables 

Most of the variables exhibit regular variations in pe¬ 
riods and amplitudes, i.e., RR Lyraes and eclipsing bi¬ 
naries. Table 2 lists the types of the 1237 variables from 
the first 2-year survey data, which are classified based 
on the shapes of their light curves, amplitudes and pe¬ 
riods of phase diagrams. Half of the detected variables 
are found to be RR Lyres according to their specific pe¬ 
riods and light curve shapes. The remained variables 
include 431 eclipsing binaries, 72 Semiregular pulsators, 
29 Mira stars, 11 slow irregular variables, 11 RS Canum 
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TABLE 2 

Distribution of Types of Variables from the TNTS 


Variable Type 

N 

% 

N(New) 

RR lyr 

661 

53.44 

59 

Eclipsing binaries 

431 

34.84 

157 

Semiregular pulsators 

72 

5.82 

68 

Mira (LPV) 

34 

2.75 

0 

Irregular 

11 

0.89 

0 

RS Canum Venaticorum 

11 

0.89 

2 

Gamma Doradus stars 

7 

0.57 

7 

W Virginis 

3 

0.24 

0 

Delta Scuti Stars 

3 

0.24 

3 

Anomalous Cepheids 

2 

0.16 

1 

Cepheids 

1 

0.08 

1 

Nova 

1 

0.08 

1 


Venaticorum stars, 7 Gamma Doradus stars, 5 long pe¬ 
riod variables, 3 W Virginis stars, 3 Delta Scuti stars, 
2 Anomalous Ceplreids, 1 Cepheid, and 1 nove-like star. 
Note that our tentative classifications of variable stars 
are based only on the unfiltered photometry. 

4.2.1. RR Lyrae Stars 

RR Lyrae stars, located at intersections of the insta¬ 
bility strip and horizontal branch, belong to a subclass 
of pulsating variable stars. They are old, low-mass, and 
metal-poor ’’Population II” stars which are usually found 
in globular clusters. According to their pulsation modes, 
the RR Lyrae stars can be divided into three subtypes: 
RRAB-type (pulsating in the fundamental mode); RRC- 
type (pulsating in first overtone); and RRD-type (pul¬ 
sating simultaneously in two modes, Smith 1995). Mean¬ 
while, different subtypes of variable stars have different 
periods and amplitudes. In our sample, there are 526 
RRAB-type stars, 126 RRC-type stars and 9 RRD-type 
stars in light of the periods and amplitudes of variations, 
which are shown in Figure 8. For our sample, the mean 
period of the RRAB-type variables is 0.57 d, while the 
RRC-type ones have the mean period of 0.34 d. The 
mean amplitude of the light variation for the RRAB sub¬ 
group (0.84 mag) is relatively larger than that of the RRC 
subgroup (0.44 mag). 

The light curves of RRAB Lyrae stars are character¬ 
ized by variations with a fast rising followed by a slow 
declining (Jurcsik et al. 2009; Cacciari et al. 2005). The 
variable star S100-3820 (see Fig. A10) was originally de¬ 
tected by AASVO VSX and identified as a RRAB-type 
star. However, our observation shows that it has a pe¬ 
riod of ~0.5 d and a variation amplitude of 0.37 mag, 
lying at the boundary between RRAB-type and RRC- 
type (see Figure 8). According to our result, its rise time 
is equal to the decline time, which likely just evolved 
off the zero-age horizontal branch (ZAHB) and exhibits 
light variation with a fixed amplitude and a longer pe¬ 
riod (Cacciari et al. 2005, §3.2). Another two interesting 
RR Lyrae stars in our sample are S037-10947 (see Fig. 
A10) and S585-9662 (see Fig. A5), which lie at a similar 
position as S100-3820 in Figure 8. Although they have 
typical features of this type of variable stars, i.e., fast ris¬ 
ing and slow declining, they show a little bump close to 
the minimum light. More data are needed to determine 
the properties of these two interesting RR Lyrae stars. 

Some RR Lyrae stars are known to have “Blazhko ef¬ 


fect” which exhibits modulations in their periods, pulsa¬ 
tion amplitudes, light curve shapes, and radial velocities 
(Blazko 1907; Buchler & Kollatli 2011). Previous surveys 
have found that 10%-30% of the RRAB-type variables 
have such a effect (Alcock et al. 2003) and this ratio can 
reach up to 50% according to a recent study by Jurcsik 
et al. (2009). Fourteen RRAB-type Lyrae stars of our 
sample are found to have phase modulations and can be 
tentatively classified as RRAB stars with the Blazhko 
effect. Obviously, we did not detect such a high ratio 
from our sample when compared with previous surveys. 
It is likely the current TNTS data may not be qualified 
enough to see most of the phase modulations. Increasing 
the cadence of the survey may help to discover more RR 
Lyrae stars with Blazhko effect. 

The period-amplitude diagram for RR Lyrae stars, also 
termed as Bailey diagram, was used to investigate their 
properties in globular clusters more than a century ago 
(Bailey 1902). Inspecting our RRAB stars shown in Fig¬ 
ure 8, we find that these variables have different locations 
in this diagram. For example the periods could differ by 
~0.1-0.2 day at given amplitude. There has been debate 
on whether Oosterhoff type or metallicity of the clus¬ 
ters accounts for these differences. The Oosterhoff-type 
classification of the Galactic globular clusters was based 
on the frequency distribution of their RR Lyrae vari¬ 
ables’ periods (Oosterhoff 1939). For example, RRAB 
stars in Oosterhoff I (Ool) clusters have mean periods 
of ^0.55 day, while those in Oosterhoff II (OoII) clus¬ 
ters have mean periods of ~0.64 day (i.e., Smith et al. 
1995, Clement et al. 2001, Catelan 2009). On the other 
hand, the metallicity has been found to correlate with 
both period and amplitude for the RR Lyrae stars (i.e., 
Marconi et al. 2015 and references therein). It has been 
found that the observed period of RR Lyrae stars tend to 
become shorter with increasing metallicity (i.e., Preston 
1959; Kunder & Clraboyer 2009; Fiorentino et al. 2015), 
while the pulsation amplitude decreases as the metallic¬ 
ity increases (i.e., Bono et al. 1996). Such an amplitude- 
period-metallicity relation can be used to estimate metal¬ 
licity for different RRAB Lyrae stars (Sandage 1982, Al¬ 
cock et al. 2000, Sandage et al. 2004, Kinemuchi et 
al. 2006). The metallicity dependence of the period- 
amplitude relation can be interpreted as an opacity ef¬ 
fect in stellar evolution (i.e., Bono et al. 1996), but a 
detailed explanation is beyond the scope of this paper. 
Note that the effect of Oosterhoff classification on the 
period-amplitude relation may be in part coupled with 
that of metallicity, as Ool clusters tend to have higher 
metallicity than clusters of OoII (Catelan 2009). 

Note that the amplitudes of RR Lyraes shown in Fig¬ 
ure 8 are measured in the unfiltered light curves, which 
are found to be systematically smaller by 0.15 mag than 
those measured in the V band (i.e., Zorotovic et al. 
2010). A similar offset was noticed in the unfiltered 
data from other surveys such as the CSS (Drake et al. 
2013) and SSS (the Siding Spring Survey; Torrealba et 
al. 2015). We thus convert the unfiltered amplitude from 
the TNTS to the V-band value by adding an offset of 0.15 
mag for each star of our sample. The corrected period- 
amplitude diagram of the RRAB-type stars is shown in 
the inset of Figure 8, where we add the reference line 
of Ool (solid line) based on eq. (11) in Zorotovic et al. 
(2010). Overplotted is the OoII line (dashed line) which 
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follows the same trend as Ool but shifted in period by 
A LogP=+ 0.06. Following Miceli et al. (2008), a period 
shift of 0.05 d relative to the Ool line was adopted to di¬ 
vide our RRAB Lyrae stars into the Ool and OoII types. 
One can see that the majority of the TNTS RRAB Lyrae 
stars cluster around the Ool line, which is consistent with 
that inferred from the CSS sample extending up to ^60 
kpc above the galactic plane (see Fig. 5 of Drake et 
al 2013) and also consistent with that from the Siding 
Spring Survey (see Fig. 11 of Torrealba et al. 2015). 

Since we have a relatively large sample of RRAB Lyrae 
stars locating from the disk to the halo of the Galaxy, 
we can carry on studies on the spatial distributions of 
their metallicities and this helps trace history of the 
Galaxy formation (Catelan 2009; Pietrukowicz 2014). 
The metallicity is estimated using the empirical period- 
amplitude-metallicity relation given by Sandage (2004), 
which follows as: 

[Fe/H] = -(1.453 ± 0.027) A v - (7.990 ± 0.00l)logP 
-(2.145 ±0.025) (7) 

Where Ay is the V-band amplitude corrected from the 
unfiltered value, P represents the period. The amplitude- 
period-metallicity relation gives an estimate of the metal¬ 
licity that is generally consistent with the spectroscopic 
method, with an uncertainty of 0.3 dex (see Fig. 9 in 
Kinemuchi et al. 2006). 



Fig. 8. — The amplitude of light variation as a function of the 
period for the TNTS RR Lyrae stars. The circles and triangles 
represent the RRAB-type stars (fundamental overtone) and RRC- 
type stars(first overtone), respectively. The blue points and red 
points represent RRD-type stars and RRAB stars with Blazhko 
effect. The inset panel shows the Bailey diagram of the RRAB- 
type stars, with the unfiltered amplitudes converted into the V- 
band values. The solid and dashed lines represent the locations of 
the Oosterhoff I and Oosterhoff II subtypes (Zorotovic et al. 2010). 


The distance to the RRAB Lyrae stars can be calcu¬ 
lated using the relation given as: 

d = l()[ ( Y 0 h- M v+5]/5 (g) 

where (Fo)s represents the static average magnitudes in 
the V band. The Galactic extinctions of the TNTS un¬ 
filtered magnitudes were corrected from the recalibrated 


infrared dust maps (Schlafly & Finkbeiner 2011). For the 
RR Lyrae sample in common, we found that the TNTS 
unfiltered magnitudes have an offset of ~0.1 mag rela¬ 
tive to the V-band magnitudes from Drake et al. (2013). 
This difference was corrected in calculating the distances 
to the TNTS RRAB Lyrae stars. Assuming an averaged 
absolute magnitude My =0.6 mag for the RRAB Lyrae 
stars (Keller et al. 2008), we can derive their heliocentric 
distances from the TNTS data. Then we converted these 
distances into those relative to the Galactic plane of the 
Milky way (i.e., |Z|). 

Fig. 9 shows the spatial distribution of the metallic¬ 
ity derived for the TNTS RRAB Lyraes. The measure¬ 
ments from the CSS photometry and the SDSS spectra 
are overplotted for comparison. As it can be seen, the 
photometric metallicities from the TNTS and CSS obser¬ 
vations are consistent with the spectroscopic estimates 
from the SDSS sample. The mean metallicity from the 
TNTS, the CSS and the SDSS is -1.59 dex, -1.57 dex, 
and -1.64 dex, respectively. The inset of Figure 9 shows 
the TNTS RRAB Lyraes distinguished by the Oosterhoff 
types. From Figures 8 and 9, one can see that the Ool 
population can occur at larger distances perpendicular 
to the Galactic plane (i.e.,|Z|~45 kpc), and they have 
higher mean metallicity (i.e., [Fe/H] ~ —1.5 dex) than 
the OoII population (i.e, [Fe/H] ~ —2.1 dex). The metal¬ 
licity of RRAB Lyraes tends to have narrower distribu¬ 
tion along the vertical distances away from the Galactic 
plane. This tendency can be explained in part with that 
the RRAB Lyraes in Ool clusters dominate at larger ver¬ 
tical distances, while a considerable fraction of these stars 
come from the OoII clusters (i.e., ~22% for the TNTS 
data) locating at relatively smaller vertical distances (i.e., 

|Z| < 20 kpc). This dichotomy in metallcity distribu- 
tion/Oosterhoff classification of RRAB Lyrae stars sug¬ 
gests that the Galactic stellar halo was perhaps formed 
by two distinct populations, favoring the dual-halo mod¬ 
els (i.e., Miceli et al. 2008, An et al. 2013). 



0 10 20 30 40 50 60 70 


IZI (kpc) 

Fig. 9. — The distribution of the metallicity measured for the 
RRAB stars versus their vertical distances from the galactic plane 
of the Milkyway. The blue, green, and grey points represent the 
TNTS, SDSS, and CSS data, respectively. The red and blue points 
in the inset panel represent the Ool and OoII subsamples of RRAB 
Lyraes from the TNTS. 















Variables from TNTS 


7 


4.2.2. Eclipsing Binary Stars 

Based on the variations of light curves, the eclipsing 
binaries can be classified into three broad categories: 
contact (EC), semidetached (ES), and detached (ED) 
(Paczynski et al. 2006). The EC binary represents a 
typical W UMa system (EW), which shows a contin¬ 
uous change in brightness and has an almost equal or 
un-obviously varying depth between the primary and 
secondary minima. The ES binary represents the Beta 
Lyrae-type system (EB), which shows a shallower change 
in brightness during the phase out of eclipse than the EC 
binary and has deeper primary depth than that of the 
secondary. The ED binary represents Algol-type system 
(EA), which shows almost constant brightness during the 
phase out of eclipse. Among our sample, there are 313 
ECs, 88 ESs, and 30 EDs. For the 299 new variables, one 
third belongs to the eclipsing binary stars. 

There are 67 interesting EC or ES binaries which ex¬ 
hibit the O’Connell effect that the two successive out-of- 
eclipse maxima are unequally high in the light curves 
(O’Connell 1951; Nataf et al. 2010). Models such as 
captured circumstellar material (Liu & Yang 2003) and 
starspot activity (Bell et al. 1990) have been proposed to 
explain the O’Connell effect, but these explanations are 
still controversial. More photometric and spectroscopic 
observations are needed to understand the origin of such 
an effect seen in the eclipse binaries. 


4.2.3. Other Types of Variables 

The long periodic variable stars include Miras, semireg¬ 
ular pulsators, and slow irregular variables. Mira vari¬ 
able shows large-amplitude variabilities (i.e., >2.5 mag 
in V band) with periods up to several hundred days. 
Semiregular pulsators show diverse periodicities in a 
range from several days to several thousand days; and 
the amplitude in light variations can vary from 0.01 mag 
to several magnitudes 10 . In our sample, 29 are Mira 
stars, 62 are semiregular pulsators, and 11 are slow ir¬ 
regular variables. 

Observationally, Gamma Doradus stars show multiple 
periods ranging from 0.3 d to 3 d with a typical variation 
amplitude of less than 0.1 mag (Cuypers et al. 2009). 
Based on our determinations of the periodicity, 5 new 
ones (S070-6584, S598-3499, S619-728, SC49-6456, and 
SE61-2890) from our discovery perhaps belong to this 
class. Delta Scuti variables are late A- or early F-type 
pulsating stars. Their pulsation periods are found to be 
in the range of 0.02-0.25 days (Breger 2000), such as star 
S585-17826 listed in Table 1. W Virginis (CWA) is type 
II Cepheid which is regarded as the low mass analogues 
of the Cepheid (Sterken & Jaschek 1996). SK39-3665 
may be such a star in our sample. 


5. CONCLUSION 

We have presented the analysis of the unfiltered CCD 
images from the Tsinghua-NAOC transient survey. A 
total of 600 sky fields (~1300 square degrees), with re¬ 
peated observations larger than 40 times during the pe¬ 
riod from October 2012 to August 2014, were selected 
for the detections of variable stars. We carried out pho¬ 
tometry using SExtractor. Fluxes of all other images 
are normalized to the reference ones; and the resulting 
unfiltered magnitudes were further calibrated using the 
R-band magnitudes of the PPMX catalogue. We used 
SRPAstro, VARTOOLS and PDM to create light curves 
and detect variable stars. 

We totally detected 1237 variables (including 299 new 
ones) with magnitude variations of > 0.1 mag and peri¬ 
ods ranging from O.lcl to 500d. These variables include 
RR Lyraes, eclipsing binaries, and semiregular stars etc.. 
Among them 10 RRAB stars show the Blazhko effect and 
64 eclipse binaries show the O’Connell effect. 

We further examined the spacial distribution of the 
photometric metallicity for the RRAB stars and found 
that it shows a larger scatter near the Galactic plane 
(i.e., —3 < [Fe/H] < 0) but tends to converge at larger 
distances (i.e., —2 < [Fe/H] < —1. This variation can 
be explained with that the Ool type of RRAB Lyrae 
stars have a higher metallicity and can extend out to a 
larger vertical distance (i.e., |Z| ~ 45 kpc) while the OoII 
type of RRAB stars have a lower metallicity and they are 
relatively abundant at a vertical distance with |Z| < 20 
kpc. This result favors that the Galactic halo was formed 
by two components with distinct properties. 
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APPENDIX 



TABLE 1 

VARIABLE STARS 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) (7) 

(8) 

(9) 

(10) 

(ID 

(12) 

(13) 

(14) 

Field 

ID 

R.A. 

Dec. 

Unfilter 

Av Amplitude 

Js 

Period(TNTS) 

Src 6 

Period(VSX) 

rpc 

1 o 

Type d 

Note 6 



(J2000.0) 0 


(mag) 


(d) 


(d) 

(HJD- 2456000) 




S007 

424 

22:47:55.16 

+25:46:59.60 

13.88 

0.19 

0.33 

0.89 

0.382994 

(LS) 

0.382989 

196.824 

EC 

V,C 

S008 

1839 

22:51:42.89 

+25:32:21.01 

14.59 

0.35 

0.45 

1.59 

0.354250 

(LS) 

0.354265 

196.650 

RRC 

V,C 

S008 

542 

22:53:00.98 

+25:43:48.03 

14.84 

0.34 

0.23 

1.23 

46.094477 

(LS) 


214.880 

SR 

N 

S008 

729 

22:53:18.25 

+25:42:40.00 

14.35 

0.35 

1.11 

8.48 

0.471269 

(LS) 

0.471272 

196.318 

RRAB 

V 

S008 

8407 

22:54:12.33 

+24:38:07.03 

15.41 

0.26 

1.02 

2.46 

0.608931 

(LS) 

0.608926 


RRAB 

v,c 

S008 

8461 
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SR 

N 

S009 

9362 

23:01:39.04 

+24:26:04.04 
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15.57 
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23:35:53.93 
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15.31 
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0.497227 
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14.44 

0.35 

1.10 

4.87 

0.471307 

(LS) 
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N 


Note:- 

a : Based on the TNTS reference images. 

b : L-S: Lomb-Scargle method; PDM: Phase Dispersion Minimization method; VSX: from VSX catalogs. 
c : Epoch when primary eclipse or minimum light occurs. 

d : ACEP= Anomalous Cepheids; RS= RS Canum Venaticorum; EC= contact binary; ES= semi-detached binary; ED= detached binary; CWA= W Virginis; RR= RR Lyrae 
variable; BL= “Blazhko effect”; CEP= Cepheids ; GDOR= Gamma Doradus star; DSCT= Delta Scuti star; SR= Semiregular pulsator; M= Mira variable; LPV= long period 
variable; L= slow irregular variable 

e : V= AAVSO VSX variable; C= CSS variable; L= LINEAR variable; N= New variable. 
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Fig. 1a.— Phase diagrams of new variables. 
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Fig. 9a.— Continued. 
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Fig. 10a. — Continued (Including two special RRAB stars which are not new). 










